Getting a kick from equal-mass binary black hole mergers 
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The final evolution of a binary black-hole system gives rise to a recoil velocity if an asymmetry is present in 
the emitted gravitational radiation. Measurements of this effect for non-spinning binaries with unequal masses 
have pointed out that kick velocities ~ 175 km/s can be reached for a mass ratio ~ 0.36. However, a larger 
recoil can be obtained for equal-mass binaries if the asymmetry is provided by the spins. Using two independent 
methods we show that the merger of such binaries yields velocities as large as ~ 440 km/s for black holes having 
unequal spins that are antialigned and parallel to the orbital angular momentum. 

PACS numbers: 04.25.Dm, 04.30.Db, 95.30.Sf, 97.60.Lf 



Introduction. Binary black hole systems are expected to 
be one of the strongest sources of gravitational waves and 
are therefore the subject of intense and careful investigations. 
With earth-based gravitational-wave detectors now working 
at design-sensitivity and a space-borne detector in its formu- 
lation phase, the need for reliable templates to be used in 
matched filtering techniques has inspired a renewed enthusi- 
asm in numerical analysis. Using numerical methods devel- 
oped recently |ll H B|] there has been an explosion of results 
(see, e.g., lefsjEBuUUUMMMMMMM)- 

These developments are important for at least three differ- 
ent reasons. First, they allow for improved templates to be 
used in the analysis of the data coming from the detectors. 
Second, they allow probes of General Relativity in regimes 
that have previously been inaccessible. Last but not least, 
they can provide, even through the solution of the Einstein 
equations in vacuum, important astrophysical information. 

Together with energy and angular momentum, gravitational 
radiation also carries away linear momentum. In the case of 
a binary system of non-spinning black holes, a physical intu- 
ition of this loss of linear momentum can be built rather easily. 
As the two bodies orbit around the common center of mass, 
each will emit radiation which is forward-beamed. Unless the 
two black holes have exactly the same mass, their motion will 
be different, with the smaller black hole moving more rapidly 
and, hence, being more efficient in beaming its emission. The 
net momentum gained over an orbit is negligible if the orbit is 
almost circular (the momentum loss in any direction is essen- 
tially balanced by an equal loss in the diametrically opposite 
direction), but it can become large when integrated over many 
orbits, leading to a recoil that is a fraction (< 10~ 2 ) of the 
speed of light during the last portion of the orbit prior to the 
merger. 

A number of PN/perturbative analyses (see, e.g., lfl7l 
[l8ll ) have provided estimates of this recoil velocity, while 
numerical-relativity simulations JH [I]]] have recently mea- 
sured it to rather high precision, predicting a maximal kick 
of 175 km/s for a binary system of nonspinning black holes 
with a mass ratio q = M1/M2 — 0.36, where Mi and M2 are 
the masses of the two black holes. Such a recoil has indeed 
quite important astrophysical consequences, since it could, 



provided it is large enough, kick the binary out of its host 
environment. Clearly, a replaced or an even missing central 
black hole would have dramatic consequences for the further 
development of the host. Determining accurately what are the 
expected escape velocities for the most typical environments 
hosting a binary black hole system is rather difficult, but the 
estimates made in refs. IU9I1 . for instance, predict that the es- 
cape velocities for dwarf galaxies and globular clusters are 
< 100 km/s, but for giant galaxies these can be ~ 1000 km/s. 

When adopting a purely geometrical viewpoint, it is ob- 
vious that a kick velocity should be expected in any binary 
system which is not perfectly symmetric. A difference in the 
masses is a simple way of producing such an asymmetry but 
surely not the only one. Indeed, even an equal-mass system 
can be made asymmetric if the two black holes have unequal 
spins. Also in this case, a simple physical intuition can be 
constructed. Consider, for simplicity an equal-mass binary in 
which only one member is spinning parallel to the orbital an- 
gular momentum. As a result of the spin-induced frame drag- 
ging, the speed of the nonspinning body will be increased and 
its radiation further beamed. Using PN theory at the 2.5 order, 
Kidder [20] has treated this spin-orbit interaction concluding 
that in the case of a circular, non-precessing orbit, the total 
kick for a binary system of arbitrary mass and spin ratio can 
be expressed as 01711 
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where 121,2 = Si^/Mf 2 are the dimensionless spins of the 
two black holes and these are aligned with the total orbital 
angular momentum, i.e., Si,2 — ai. 2 M 2 2 e z for an orbital 
motion in the (x, y) plane. Here, c\ and c 2 are factors de- 
pending on the total mass of the system and on the orbital 
separation at which the system stops radiating. This radius 
is difficult to determine precisely as it lies in a region where 
the PN approximation is not very accurate and is, in practice, 
not even a constant but, rather, depends on both the mass and 
the spin ratio. Assuming for simplicity c\ ~ C2, expression 
(fl]l reveals that a substantial contribution to the recoil velocity 
comes from the spins alone. In addition, for any given q, it 



±x/M ±p/M mi/M m 2 /M Si/M 2 S 2 /M 2 M ADM /M 

rO 3.0205 0.1366 0.4011 0.4009 -0.1460 0.1460 0.9856 

rl 3.1264 0.1319 0.4380 0.4016 -0.1095 0.1460 0.9855 

r2 3.2198 0.1281 0.4615 0.4022 -0.0730 0.1460 0.9856 

r3 3.3190 0.1243 0.4749 0.4028 -0.0365 0.1460 0.9857 

r4 3.4100 0.1210 0.4796 0.4034 -0.0000 0.1460 0.9859 



TABLE I: The puncture initial data parameters defining the binaries: 
location (±x/M), linear momenta (±p/M), masses (mi/M), spins 
(Si/M) and ADM mass measured at infinity (A/ ADM ). 



predicts a linear growth of the recoil velocity with increasing 
difference in spins, yielding a kick which is comparable with 
the one coming from the asymmetry in the mass. Stated differ- 
ently, when it comes to recoil velocities, the spin contributions 
may be the dominant ones. 

Techniques. The numerical evolutions have been carried out 
using a conformal-traceless formulation of the Einstein equa- 
tions as described in [21], with "1+log" slicing and T-driver 
shift, and advection terms applied to the gauge conditions 
as suggested in J3]. Spatial differentiation is performed via 
straightforward finite differencing using fourth-order stencils. 
Individual apparent horizons are located every few timesteps 
during the evolution [22]. Vertex-centered adeptive mesh- 
refinement (AMR) is employed using nested grids l23ll with 
the highest resolution concentrated in the neighbourhood of 
the individual horizons. For each of the models studied, 
we have carried out simulations with fine-grid resolutions of 
h = 0.030 M and h = 0.024 M, where M = M x + M 2 is 
the total mass of the system. For a subset we have carried out 
a further evolution at h = 0.018 M to assess the convergence 
and estimate the error for the h — 0.024 M results. 

The initial data is constructed using the "puncture" 
method ll24ll . which uses Bowen-York extrinsic curvature and 
solves the Hamiltonian constraint equation numerically as 
in ||25jl. We have considered a sequence of binaries for 
which the initial spin of one of the black holes is fixed at 
S2/M 2 — 0.146e z , and, to maximize the recoil [cf. eq. (Q]l], 
the spin of the second black hole is in the opposite direction 
and has a modulus which is varied in steps of 1/4. As a re- 
sult, we obtain 5 sets of binary black holes having spin ratios 
— 1, —3/4, . . . , 0. The orbital parameters (m», X\ andp^) for 
spinning binaries in quasi-circular orbit are determined using 
an effective potential method ll26tl (a minimum is found in the 
binding energy of the system) in such a way that the black hole 
masses (Mj = y/M? + Sf/{AM?)) are equal and M = 1. 
The data for the initial separation and momenta are summa- 
rized in Table U and are chosen so that all of the binaries have 
the same orbital angular momentum. 

To avoid possible systematic errors and improve the accu- 
racy of the measurements, the kick velocity has been com- 
puted using two different and independent methods. The first 
and more traditional one makes use of the Newman-Penrose 
quantity ip4 to provide the rate of change of the linear momen- 
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FIG. 1: Recoil velocity as function of time for a binary system of 
nonspinning black holes with a mass ratio of 2/3 at an initial sepa- 
ration 4.1 M . The set of curves (a) and (b) differ in the choice of the 
integration constant, while the solid and dashed lines show the two 
independent computations of the momentum flux [eqs. (O and (3J|]. 

turn in the i-direction as i27l, l28tl 




The second and novel method, instead, uses a perturbative 
wave-extraction procedure 12911 that calculates gravitational 
waves in terms of gauge-invariant, even Q&1 and odd-parity 

Qfyl metric perturbations of the Schwarzschild spacetime af- 
ter they have been decomposed into spherical harmonics of 
spin weight —2 [f30h . Using these quantities it is possible to 
reconstruct the momentum fluxes in the x and y-directions as 

00 

P x + iP y = Y: F(Qt°\Qt o) *> Qt°\Qt o) *) , (3) 

where the overdot indicates a (coordinate) time derivative and 
the * complex conjugation [the complete expression for the 
formula ([3]l will be presented in a longer paper 13111 1. 

We have validated both methods by measuring the recoil 
velocity for a binary system of nonspinning black holes hav- 
ing a mass ratio of 2/3 at an initial separation of 4.1 M. The 
results of this calibration extracted at r = 50 M are shown in 
Fig. Q] which reports the evolution of the kick velocity using 
*5>4 (solid lines), and the gauge-invariant quantities when the 
summation in (O is truncated to the first 6 multipoles (dashed 
lines), which we have found to be sufficient to reach conver- 
gence. Indicated with symbols are the estimates and relative 
error bars obtained by |8|] (circle) and by 111 ill (star). 

We note that because the binary system starts evolving at 
a finite separation, it will have already gained a net linear 
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FIG. 2: Recoil velocity as function of time for the sequence of runs 
i.e., from rO with — ai = a2 = 0.586, to r4 with a\ = 0, a2 = 
0.586). Note that the merger is delayed for smaller values of \a\ \. 
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FIG. 3: Kick velocities and error bars for different spin ratios; the 
dashed lines show a linear fit of all the data when the point at 
a\/a,2 = 1 is given an infinite weight since [u|kick = for ai = a2. 

found to be the smallest needed distance for consistent results. 



momentum which can influence the value of the final kick. 
Computing this initial linear momentum amounts to selecting 
a proper constant in the integration of (O or (01. Fortunately, 
this is rather straightforward to do and amounts to determining 
the direction in 3-space in which the center of mass of the sys- 
tem is moving initially. In practice, we plot the evolution of 
the x and y-components of the kick velocity (the z -component 
is zero because of symmetry) and calculate the vector to the 
centre of the spiral generated as the evolution proceeds. This 
vector is then composed with the final one, yielding the final 
kick; note that being a vector this integration constant is not 
simply an additive constant for the kick velocity |w|kick- In 
Fig. [T| we have plotted the effect of including this constant, 
comparing the case where it is set to zero [set of curves (a)], 
with a value set by extrapolating the recoil backwards to com- 
pensate for the small but nonzero initial linear momentum [set 
of curves (b)]. In the first case we find agreement with IfXlTI - 
while in the second case the good agreement is with f^. 

A validation of this procedure is also rather straightforward: 
only an accurate estimate of the initial momentum yields a 
monotonic evolution of the kick velocity (or, in the case of 
very close binaries, reduces the oscillations considerably); any 
different choice would yield the oscillations seen in curves (a) 
(cf, Fig. 1 of ref. H or Fig. 3 of ref. HI). Clearly, select- 
ing the correct integration constant becomes less important as 
the separation in the binary is increased (see also the discus- 
sion below), but it can easily lead to errors of 10% or more 
for the rather close binaries considered here. A more detailed 
discussion on the integration constant will be given in 113 ill . 
Results. In Fig.|2]we show the evolution of the recoil velocity 
for the 5 binaries considered, with the dashed lines referring 
to the gauge-invariant extraction and the solid lines to the one 
using ^4; both quantities are extracted at 50 M which we have 



A number of interesting features are worth remarking. 
Firstly, all of the curves show a monotonic growth as a result 
of a suitable choice for the integration constant. Secondly, 
the two types of measurements agree to very good precision. 
Thirdly, binaries that have a spin ratio closer to zero merge 
progressively later. In the case of the binaries considered here 
this effect is partially masked by the fact that the binary rO is 
at a smaller separation than r4 (cf., Table H}. Bearing this in 
mind, it is however apparent that the growth rate of the kick 
velocity (and hence the rapidity of the inspiral), increases with 
the asymmetry in the spins. Fourthly, increasing the initial 
separation for a binary with ai/a-2 = —1 does not change 
significantly the integration constant chosen for rO, thus indi- 
cating that the kick estimate for the latter is robust. Finally, as 
in unequal-mass binaries, the largest contribution to the kick 
comes from the final parts of the inspiral and is dominated by 
the last orbit. However, unlike equal-mass binaries, the post- 
merger evolution of the kick velocity is not modified substan- 
tially by the quasi-normal mode ringing (cf., Figs. [2] and Q]), 
with the final kick velocity being only slightly smaller than 
the maximum one reached during the evolution. 

As predicted by the PN expression (Q~|), the final velocities 
shown in Fig.|2]exhibit a linear dependence with the spin ratio, 
and this is shown in Fig. [3] which reports the asymptotic kick 
velocities when measured with '5 4 (open circles) or with the 
gauge-invariant perturbations (stars). Also indicated are the 
error bars which include errors from the determination of the 
integration constants, from the dependence of the waveforms 
on the extraction radii, and from the truncation error. 

The data points in Fig. [3] are not the only ones available and 
indeed a binary system with ai /a,2 — 1 is bound to produce 
a zero kick velocity. The dashed lines in Fig. [3] represent a 
linear fit of all the data when the point at a\ ja% — 1 is given 
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an infinite weight to account that |f|kick = when a\ = 
(short-dashed line for ^4 and long-dashed for Q^^)- These 
lines are only indicative and bear a physical significance only 
if the linear dependence should hold for all the possible values 
of the spin-ratio. 

Conclusions We have calculated the recoil velocity from the 
inspiral and merger of binary black-hole systems with equal 
masses but unequal spins. To increase the accuracy and re- 
move systematic errors, we have performed our measure- 
ments using both the Newman-Penrose quantity ^4, as well 
as gauge-invariant metric perturbations of a Schwarzschild 
spacetime. The two methods agree to very good precision and 
indicate that the recoil velocity produced by spin-asymmetries 
can be considerably larger than that expected from asymme- 
tries in the masses. More specifically, we have found a max- 
imum recoil velocity of 257 ± 15 km/s for a system having 
a spin ratio a\ja2 = — 1 and = 0.584. In addition, a 



linear scaling holds between the kick velocity and the spin 
ratio, thus qualitatively confirming the PN expectation. A 
linear extrapolation to the case of extremal black holes, i.e., 
a\ = 1 = — a-2, suggests a maximal value for the kick veloc- 
ity from unequal-spin, but equal-mass binaries of ~ 440 km/s. 
Such a recoil velocity would be more than twice that produced 
by non-spinning but unequal-mass binaries. 

At the time of the submission of this paper a number of 
complementary studies have also been submitted I132ll33ll34ll . 
which support the results found here as well as examine non- 
aligned spin configurations. 
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